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CAS A: THE BRIGHT X-RAY KNOTS AND OXYGEN EMISSION 

D. Dewey, 1 T. DeLaney, 1 and J.S. Lazendic 2 
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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. Chandra High-Energy Transmission Grating (HETG) X-ray spectra are extracted from 
17 bright, narrow regions of CasA and provide unique measurements of their kinematic and plasma states. 
From the dominant emission lines, e.g. He-like Si, we derive accurate Doppler shifts in the range —2500 to 
+4000 kms -1 ; these agree well with transverse- velocity measurements and allow the features to be located in 
3D. Plasma diagnostics of these regions indicate temperatures largely around IkeV with some above 3keV. 
Using as well the non-dispersed zeroth-order data, we determine NEI model parameters for the regions which 
lead to density estimates. Values of n e « 100 cm -3 are likely the maximum of a range of densities in this X- 
ray emitting material. The common "oxygen-rich" assumption is coarsely tested by comparing the integrated 
OVIII line flux and continuum levels. It appears that most of the continuum is due to another source, e.g., 
from He. 

ABSTRACT 

Chandra High-Energy Transmission Grating (HETG) X-ray spectra are extracted from 17 bright, narrow 
regions of CasA and provide unique measurements of their kinematic and plasma states. From the dominant 
emission lines, e.g. He-like Si, we derive accurate Doppler shifts in the range -2500 to +4000 kms -1 ; these agree 
well with transverse- velocity measurements and allow the features to be located in 3D. Plasma diagnostics of 
these regions indicate temperatures largely around 1 keV with some above 3 keV. Using as well the non-dispersed 
zeroth-order data, we determine NEI model parameters for the regions which lead to density estimates. Values 
of n e « 100 cm -3 are likely the maximum of a range of densities in this X-ray emitting material. The common 
"oxygen-rich" assumption is coarsely tested by comparing the integrated O VIII line flux and continuum levels. 
It appears that most of the continuum is due to another source, e.g., from He. 

Key Words: RADIATION MECHANISMS: THERMAL — SUPERNOVA REMNANTS: INDIVID- 
UAL(CAS A (Gill. 7-2.1)) — TECHNIQUES: SPECTROSCOPIC — X-RAYS: ISM 



1. HETG OBSERVATION OF CASA 

The X-ray emitting material in young super- 
nova remnants is the bulk of the ejected mass and 
gives direct information on the progenitor, the ex- 
plosion mechanism, and illuminates the overall hy- 
drodynamics at the present epoch. Cas A is a bright, 
Galactic remnant which allows detailed study in all 
wavebands. Here we add high spatial and spectral 
resolution X-ray measurements to the CasA knowl- 
edge base. 

The HETG is a slitle ss, dispersive spectrome- 
ter (|Canizares et al.l 12005) — nominally designed for 
point source observations. In contrast, CasA is a 
very extended source compared to the HETG disper- 
sion scale, Figure [TJ the spacing between the bright 
Si dispersed images is not much larger than Cas A it- 
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self. Note that streaks of intensity variation are seen 
along the dispersed images, the most prominent low- 
intensity streaks form a wide "X" patterrjf) and are 
likely due to a region of high local absorption (indi- 
cated on Fig. [3) 

In spite of Cas A's large size, we can obtain infor- 
mation from these disp ersed d a ta usi ng general tech- 
niques as described in iDewevI ( 20021 ). In particular, 
detailed analyses were carried out fo r Cas A's bright 
Si line emission and are presented in Lazendic et al.l 
fl2006h . hereafter Paper I. In this present work we 
summarize those results in §[2]& §[3j and discuss fur- 
ther aspects in §[4] & §0 Finally, future plans in 
these areas are discussed in §[6l 

2. KNOT VELOCITIES 

To the HETG, CasA appears as a multitude of 
sources; for analysis, we selected 17 small, bright 



3 These streaks were brought to our attention by Peter Ford 
of the ACIS science team — he initially feared that they may 
be instrumental in origin! 



1 



DEWEY, DELANEY, & LAZENDIC 




Fig. 1. Dispersed image of Cas A on the ACIS-S array. The undispersed "zeroth-order" image of Cas A is seen near the 
center; with axes roughly N to left and E down. Dispersed images of Cas A in the bright Si emission are visible to either 
side of the central image. The MEG spectrum extends from upper-left to lower-right and the HEG from lower-left to 
upper-right. Five vertical low- intensity regions indicate the gaps between the six ACIS CCD chips. The red-to-blue 
color range corresponds to 0.5 keV to 2.3 keV. 
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Fig. 2. The 17 regions are indicat ed on an X-ray i mage of 
Cas A (J2000 coordinates), from lLazendic etaD ([20 



3. PLASMA PROPERTIES 

The f jr and H- like /He- like line ratios were used 
with an NEI-model-generated line-ratio grid to de- 
termine values of kT and ionization age, r = n e £, 
for the regions studied. These values can be dis- 
played on a "Laming plot" , Figur e [H Note that un- 
like th e knot series measured in Laming fc Hwangl 
(|2003l ). the knots here were selected for their bright- 
ness and small size and so we don't expect them to 
trace out the full extent of a model curve but rather 
they should be located at the places of maximum 
Si-line emission along the model curves. The three 
high-temperature regions, R8, R10, and R17, clearly 
stand out from the general model expectations; R17 
appears to be involved with the NE jet phenomenon, 
while R8 and R10 do not at first glance show a reason 
for their high temperatures. 

B 



knots (here, "knots" is used interchangeably with 
other terms: regions, features, filaments...), as in- 
dicated on Figure O We used the filament- analysis 
technique (Paper I) to narrow the features and in- 
crease the spectral resolution. The custom- created 
PHA and RMF files were then fit in ISIS (jHouckl 
2002f ) to obtain line wavelengths and flux ratios. 



Analysis of the line shifts of the bright He-like Si 
lines gives accurate Doppler velocitie s, in agreement 
with t he transverse proper motions of DeLanev et al.1 
(|2004l ). Assuming a uniform-expansion velocity 
field, the regions can be assigned a line-of-sight z- 



coordinate: z = const x V m 



and the regions 



can be located in 3D, as shown in Figure 03 The 
trend to blue-shifted velocities in the SE agrees qual- 
itatively with early dispersive spectroscopy results 
([Markert et al.lll983h . 
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Fig. 3. The regions in 3D. Bottom: cartoon version of the 
regions on the sky. Top: the same regions viewed from 
"above" with the axes E-W and Front-Back indicated. 



REVMEXAA(SC) CSM AND LATE STAGES. 



3 



10 s 



Hn o 



7^ 



> R15 



R 1 3 



10 9 1Q 10 1Q 11 1Q 12 

onlzation age (s/cm A 3) 



10 13 



Fig. 4. Regions located in T vs. r space; dotted lines 
indicate the min-max rang e of th e model expectations 
shown in Lam ing fc Hwanj (l2003h . 



Fixing the kT and r values, we then fit the 
zeroth-order spectrum of each region in the range 
above 1.1 keV to determines the model "norm" 
and abundance values. The continuum in the 
model fitting was assu med to come from oxygen 
( Laming fc Hwangll2003h , see however the discussion 
in §0 Using a distance to Cas A of 3.4 kpc and esti- 
mated volumes for the regions, the NEI model values 
can be converted to densities of the ions and elec- 
trons in the regions (Paper I, Appendix B.) 

The derived properties of these regions are plot- 
ted as a function of their 3D radius from the ex- 
pansion center, Figure [SJ There does not appear to 
be any clear trend of temperature or density with 
this radius as simple models (next section) might 
suggest. The two outer-most regions, R16 and R17, 
show solar-like abundances of Mg whereas the inner 
regions have Mg at generally sub-solar levels. In the 
plot of time-since-shocked, t s hock = T/n e , the back- 
side regions closer to the center show lower values 
suggestive of their having more recently encountered 
the reverse shock (RS). 

4. EJECTA DENSITY DISTRIBUTION 

The eject a densities of the measured regions are 
of order 100 cm -3 ; assuming this level is typical 
throughout the full 3D shell from 100" to 130" radius 
leads to a very large mass estimate of the continuum- 
providing oxygen (Paper I.) How do these inferred 
densities compare with expected densities and the 
generally assumed model picture? 
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Fig. 5. Properties of the regions plotted vs. r3D- Unless 
otherwise indicated, the "*" and diamond symbols cor- 
respond to regions on the back and front of the remnant, 
respectively. Lines on the plots serve to "guide the eye." 
(For clarity, the outlying region R17 is plotted at a radius 
of 152" instead of its actual value of 172" .) 

The high temperature X-ray gas should trace out 
the bulk hydrodynamics of the remnant, as opposed 
to the generally un-coupled optical FMKs and the 
more mysteriously coupled (magnetic fields and rel- 
ativists elect rons) radio emission. Th e self-similar 
solution, e.g. J Chevalier fc Oishil (2003), for an s = 2 
"wind" circumstellar medium shows a range of densi- 
ties in the shocked eject a with highest density at the 
contact discontinuity (CD, see their Figure 2.) Such 
a trend of density with radius is not seen in our data. 
However, as suggested by Laming (private commu- 
nication), the radial location of the RS and CD may 
vary around the remnant and so the effective "thick- 
ness" of the high-density shell is much smaller than 
simply the range of radii over which we see high 
densities. Further, it is likely that Ray leigh- Taylor 
(R-T) instabilities will distort the otherwise sharp 
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CD ([Chevalier et all Il992l : fowarkadasl l2000h : note 
that even with R-T effects a similarly large range of 
eject a densities may still be maintained, Figure [U 
The bright, dense regions measured here then reflect 
the higher-density portions of the hydro structures 
developed by the ejecta between the RS and CD. 
Extrapolating from their density values to an over- 
all mass estimate requires further assumptions of the 
underlying density distribution. 




Fig. 6. Histograms of model ejecta densities. The images 
at upper left and right show cross-sections of a simple 
self-similar distribution and an R-T-hydro density dis- 
tribution, respectively. Histograms of the densities in a 
3D shell from these models are plotted below; the R-T- 
hydro simulation is the solid line. Both models show a 
range of densities extending to 10 or more times the most 
common value. 

5. OXYGEN-RICH ? 

It has been noted that bright eject a- dominated 
regions of CasA cannot be modeled with a pure 
high-Z metals plasma due to a need for addi- 
tio nal continuum flux from elements with "Z < 
8" ( Hughes et alJl2QQol ). The continuum could arise 
fro m mixed-in outerlayer ma teria l: He, N, perhaps 
C (Perez-Rendon et al. 2002 ). In Laming fc Hwangl 
(|2QQ3h and our Paper I this continuum was assumed 
to come from oxygen for modeling purposes. 

One reason that the actual amount of oxygen is 
an unsettled issue is that Cas A is strongly absorbed 
with an TVh of order of 10 22 cm -2 , greatly reduc- 
ing the OVIII emission. Even so, the ACIS imag- 
ing data does include oxygen-band emission which 
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Fig. 7. Oxygen in CasA. Bottom image: the oxygen 
band, 0.5 keV to 0.8 keV, from the non-grating observa- 
tion, obsid 5196. The main O-bright regions are labeled. 
Upper image: the MEG-dispersed events on SI (binned 
and smoothed) from obsid 1046; discrete O VIII line emis- 
sion is just visible from portions of the NE, N, and NW 
regions. Vertical lines indicate the extraction region for 
the spectrum of Figure [8] 

shows greatest intensity in four general regions: la- 
beled NE, N, NW and SE in the bottom image of 
Figure 

It is possible to look for and set limits on the 
oxygen emission using the OVIII Ly a and Ly (3 
lines seen in the MEG spectrum of CasA; these 
lines were also measured with the XMM- Newton 
RGS (jBleeker et al.ll200ll ). In that work the La-to- 



L/3 ratio was measured versus cross-dispersion loca- 
tion, separating the SE, NE, N regions. We defer 
a similar detailed analysis and instead consider here 
the global spectrum from the combined N & NW re- 
gion. 
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Fig. 8. Oxygen seen with the MEG. The histogram 
extracted from the extended N fe NW region is well- fit 
by a set of XSPEC vnei models which include He, Fe, 
and O. The underlying He and He+Fe components are 
also shown; the model is for the kT = 0.7 keV and 
r = 0.2 x 10 11 scm" 3 case of Tabled 

The extracted counts spectrum using our fila- 
ment technique is shown in Figure El This spec- 
trum is fit by an absorbed vnei model containing 
only He, Fe, and O. The observed lines are broader 
than the spatial image and velocity components of 
+ 1500kms- 1 and -2100 kms" 1 are included. The 
difficulty of working in this heavily absorbed range 
with relatively low resolution data is demonstrated 
by the fitting results given in Table [U The table 
shows that comparably good fits can be obtained 
over a large (but relevant) range of kT and r, hence, 
the data really constrain only four of the six model 
parameters. Note that Nn is considered a free pa- 
rameter here both for illustration of its effect on the 
fit parameters and to allow the possibility that it dif- 
fers from the expecte d range of 1.0-1.2 x 10 22 cm -2 



([Keohane et al.lll996h . 



Even with this degeneracy some preliminary con- 
clusions are possible: i) the integrated Fe/O abun- 
dance ratio is greater than 0.2 solar and ii) there 
is substantial He (and/or C, N) compared to O-fFe, 
with M(0+Fe)< 0.01 M(He). We can conclude that 
CasA is dominated by O burning products (Si, S, 
etc.) rather than by oxygen per se. 

As a caveat, note that given the large region 
extracted and the unavoidable inclusion of a large 
radial range, which includes the blastwave, it is 
possible that a finite contribution to the mea- 
sured ^conyirmmi comes from synchrotron radia- 
tion ([Hughes et al.ll2000h . 



CasA and given some information on the oxygen 
present in the remnant. Some additional work in 
these directions is listed here: 

• It should be possible to measure the Doppler 
shifts of additional knots in this data set as well 
as study individual bright oxygen features. 

• We can compare the (sparse) 3D X-ray struc- 
ture here with t h e opt ical /IR "ring" struc- 



tures ( Reed et ah! ( 1995f ) and unpublished ad- 
ditional work.) 

• We have not yet examined the dispersed Fe-K 
emission in this data to see if it contains useful 
information; if the Fe is indeed in fuzzy bub- 
bles then our slitless dispersive data will have 
reduced spectral resolving power. 

• The non-dispersed Chandra 1 Ms data set is a 
unique resource and these grating results can 
be used to confirm/calibrate the CCD-measured 
Doppler shifts. 

• It will be useful to create a "CasA knot cata- 
log" to allow comparisons across analyses, wave- 
lengths, and epochs. 

• Although Suzaku lacks the spatial resolution to 
isolate individual knots, it does have better en- 
ergy resolution in the oxygen band and may help 
constrain global Cas A emission models. 

• X-ray features in Cas A are changing in 

time ( Patnaude fe Fesenl l2QQ6h and a future 
HETG re-observation (At > 6 years) could give 
useful information on knot evolution. 

In the modeling domain, the analyses presented 
here suggest a larger goal of building an approxi- 
mate 3D model of the CasA remnant (at a given 
epoch) from a set of physically-consistent emission 
components. Work with this flavor is already under- 
way in astrop hysics, e.g., in the modeling of plane- 
tary nebulae ( Morisset et al 112005: Steff en fe Lopezl 



20061 ). and will be accelerated by making more use 
of the 3D capabilities of current computer hardware 
and software. 

Finally, in the more distant hardware future, one 
might imagine a "soft X-ray spectrometer" with a 
long-slit configuration operating in this low-energy 
rangjfl to better isolate and measure the emission 
from complex extended sources like Cas A. 



6. FUTURE WORK 

The current HETG data has produce a small but 
tantalizing set of measurements on bright knots in 



4 A minimum range of 14.5 A to 22.5 A would include im- 
portant Fe-L and O lines. Extending shortward to include 
NeX (12 A) and/or longward to include C and N lines (26-30- 
35-42 A), would of course be desirable. 
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TABLE 1 

MODEL a FITS TO THE OXYGEN AND FE-L EMISSION FROM THE N&NW REGION 



Fixed Values Fit Parameters Derived Values 



kT 


r 




Fe/O 


O/He 


Y b 
A norm 


x 2 


n e 




r/n e 


(keV) 


(10 n scm- 3 ) 


(10 22 cm- 2 ) 






(io- 3 ) 




(cm -3 ) 


(M solar ) 


(years) 


0.7 


0.16 


1.05 


0.38 


0.017 


18.1 


1.18 


73. 


11.2 


6.9 


0.7 


0.20 


1.0 


0.46 


0.016 


14.5 


1.11 


65. 


10.0 


9.7 


0.7 


0.30 


0.9 


0.68 


0.013 


9.9 


1.08 


54. 


8.3 


17. 


0.7 


0.60 


0.8 


0.82 


0.013 


7.0 


1.10 


45. 


7.0 


42. 


0.7 


1.30 


0.7 


0.53 


0.031 


4.8 


1.08 


38. 


5.8 


109. 


0.7 


3.00 


0.7 


0.15 


0.154 


4.3 


1.05 


36. 


5.4 


266. 



0.9 


0.20 


0.9 


0.69 


0.012 


8.6 


1.08 


50. 


7.7 


12.5 


0.9 


0.60 


0.6 


1.98 


0.014 


2.2 


1.04 


25. 


3.9 


74. 


0.9 


3.00 


0.4 


0.64 


0.234 


0.78 


0.93 


15. 


2.3 


623. 


1.2 


0.16 


0.8 


1.03 


0.010 


5.0 


1.07 


38. 


5.9 


13.1 


1.2 


0.60 


0.4 


4.87 


0.016 


0.74 


0.94 


15. 


2.3 


128. 


1.2 


3.00 


0.6 


0.69 


0.290 


2.9 


1.04 


29. 


4.4 


323. 



a The model is an XSPEC wabs*vnei with kT and r fixed as tabulated. The elemental abundances are: H=l, He=1000, 
O and Fe are free, and all others are set to 0.0. 



b Anorm is the usual XSPEC model "norm" value; see Paper I Appendix B for details. 

c These physical parameters are based on an emitting volume equal to the cross- dispersion width (128") times one 
quadrant of the annulus from 100"to 130": V = 0.25 * 128 * tt(130 2 - 100 2 ) « 9.3xl0 55 cm 3 . The equations of Paper I 
Appendix B are used with a distance to Cas A of 3.4 kpc. 



In all this we are guided by our current theories, 
but the data have the last word and there may well 
be some "non-spherical-elephant" surprises along the 
way. 
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